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Existence of a weak solution to a fluid-elastic 

STRUCTURE INTERACTION PROBLEM WITH THE NAVIER 
SLIP BOUNDARY CONDITION 
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Abstract 

We study a nonlinear, moving boundary fluid-structure interaction (FSI) problem between 
an incompressible, viscous Newtonian fluid, modeled by the 2D Navier-Stokes equations, and an 
elastic structure modeled by the shell or plate equations. The fluid and structure are coupled 
via the Navier slip boundary condition and balance of contact forces at the fluid-structure inter¬ 
face. The slip boundary condition might be more realistic than the classical no-slip boundary 
condition in situations, e.g., when the structure is “rough”, and in modeling FSI dynamics near, 
or at a contact. Cardiovascular tissue and cell-seeded tissue constructs, which consist of grooves 
in tissue scaffolds that are lined with cells, are examples of “rough” elastic interfaces interacting 
with an incompressible, viscous fluid. The problem of heart valve closure is an example of a FSI 
problem with a contact involving elastic interfaces. We prove the existence of a weak solution 
to this class of problems by designing a constructive proof based on the time discretization 
via operator splitting. This is the first existence result for fluid-structure interaction problems 
involving elastic structures satisfying the Navier slip boundary condition. 


1 Introduction 

We study a nonlinear, moving boundary fluid-structure interaction (FSI) problem between a viscous, 
incompressible Newtonian fluid, modeled by the Navier-Stokes equations, and an elastic shell or 
plate. The fluid and structure are coupled through two coupling conditions: the Navier slip boundary 
condition and the continuity of contact forces at the fluid-structure interface. The Navier slip 
boundary condition states that the difference, i.e., the slip between the tangential components of the 
fluid and structure velocities is proportional to the tangential component of the fluid normal stress 
evaluated at the fluid-structure interface, while the normal components of the fluid and structure 
velocities are continuous. The main motivation for using the Navier slip boundary condition comes 
from fluid-structure interaction (FSI) problems involving elastic structures with “rough” boundaries, 
and from studying FSI problems near a contact. 

FSI problems involving elastic structures with rough boundaries appear, for example when study¬ 
ing FSI between blood flow and cardiovascular tissue, whether natural or bio-artificial, which is lined 
with cells that are in direct contact with blood flow. Bio-artificial vascular tissue constructs (i.e., 
vascular grafts) often times involve cells seeded on tissue scaffolds with grooved microstructure, 
which interacts with blood flow. See Figure [l] To filter out the small scales of the rough fluid 
domain boundary, effective boundary conditions based on the Navier slip condition have been used 
in various applications, see e.g., a review paper by Mikelic [41], and mmm- Instead of using the 
no-slip condition at the groove-scale on the rough boundary, the Navier slip condition is applied at 
the smooth boundary instead. 

Another motivation for using the Navier slip boundary condition comes from studying problems 
near or at a contact (or collision). It has been shown recently that the no-slip condition is not a 
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Figure 1: Two examples of fluid-elastic structure interaction problems where the Navier slip condition may 
be more appropriate. Left: An image of a nano-patterned film surface for vascular tissue engineering |58] 
giving rise to a “rough” surface interacting with blood flow. The grooves in this tissue construct are to be 
lined with cells, producing a bio-artificial vascular graft for a replacement of diseased arteries. The grooves 
are perpendicular to the direction of blood flow. Right: A sketch of a closed mitral valve interacting with 
blood flow. Modeling the contact between leaflets with the no-slip condition gives rise to leaky valves due 
to the no-slip condition paradox in collision of smooth bodies |'49l 1261 127', 52:. 


realistic physical condition to model contact between smooth rigid bodies immersed in an incom¬ 
pressible fluid [491126, 123 [52] . It was shown that two smooth rigid bodies cannot touch each other if 
the no-slip boundary condition is considered. One solution to this no-collision paradox is to consider 
bodies with “non-smooth” boundaries, in which case collisions can occur [T8]. The other explanation 
for the no-collision paradox is that the no-slip boundary condition does not describe near-contact 
dynamics well, and a new model and/or a different boundary condition, such as for example, the 
Navier slip boundary condition, need to be employed to model contact between bodies/structures 
interacting while immersed in an incompressible, viscous fluid |47j . Examples include applications 
in cardiovascular sciences, for example, modeling the closure of heart valves. It is well known that 
numerical simulation of heart valve closure suffers from the “numerical” leakage of blood through 
a “closed” heart valve whenever the no-slip boundary condition is used. Different kinds of “gap” 
boundary conditions have been used to get around this difficulty, see e.g., m- Considering the 
Navier slip boundary condition near or at the closure would provide a more realistic modeling of the 
problem. 

From the mathematical analysis point of view, the fist step in the direction of studying the 
Navier slip boundary condition near or at a contact was made by Neustupa and Penel m who 
proved that when the no-slip boundary condition is replaced with the slip boundary condition, 
collision can occur for a prescribed movement of rigid bodies. Recently, Gerard-Varet and Hillairet 
considered an FSI problem involving a movement of a rigid solid immersed in an incompressible 
Navier-Stokes flow with the slip boundary condition and proved the existence of a weak solution up 
to collision m- In a subsequent work |20| they proved that prescribing the slip boundary condition 
on both the rigid body boundary and the boundary of the domain, allows collision of the rigid body 
with the boundary. The existence of a global weak solution which permits collision of a “smooth” 
rigid body with a “smooth” fluid domain boundary was proved in m- For completeness, we also 
mention several recent works where the slip boundary condition was considered in various existence 
results for FSI problems involving rigid bodies and Newtonian fluids BEIHSIEZI. All the above- 
mentioned works consider FSI between rigid bodies and an incompressible, viscous fluid. To the 
best of our knowledge there are no existence results for non-linear moving boundary FSI problems 
involving elastic structures satisfying the Navier slip boundary condition. The present work is the 
first existence result involving the Navier slip boundary condition for a fluid-structure interaction 
problem with elastic structures. 

Classical FSI problems with the no-slip boundary condition have been extensively studied from 
both the analytical and numerical point of view (see e.g. o H7i nn and the references within). 
Earlier works have focused on problems in which the coupling between the fluid and structure was 
calculated at a fixed fluid domain boundary, see [16] . and Mug, where an additional nonlinear 
coupling term was added and calculated at a fixed fluid interface. A study of well-posedness for FSI 
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problems between an incompressible, viscous fluid and an elastic/viscoelastic structure satisfying 
the no-slip boundary condition, with the coupling evaluated at a moving interface, started with 
the result of daVeiga [4], where existence of a strong solution was obtained locally in time for an 
interaction between a 2D fluid and a ID viscoelastic string, assuming periodic boundary conditions. 
This result was extended by Lequeurre in E3EH3, where the existence of a unique, local in time, 
strong solution for any data, and the existence of a global strong solution for small data, was proved 
in the case when the structure was modeled as a clamped viscoelastic beam. 

D. Coutand and S. Shkoller proved existence, locally in time, of a unique, regular solution for 
an interaction between a viscous, incompressible fluid in 3 D and a 3 D structure, immersed in the 
fluid, where the structure was modeled by the equations of linear elasticity satisfying no-slip at 
the interface [13] . In the case when the structure (solid) is modeled by a linear wave equation, I. 
Kukavica et al. proved the existence, locally in time, of a strong solution, assuming lower regularity 
for the initial data [32]ES 2D • A similar result for compressible flows can be found in 03]. In m 
Raymod et al. considered a FSI problem between a linear elastic solid immersed in an incompressible 
viscous fluid, and proved the existence and uniqueness of a strong solution. All the above mentioned 
existence results for strong solutions are local in time. Recently, in [3D] a global existence result for 
small data was obtained for a similar moving boundary FSI problem but with additional interface 
and structure damping terms. 

In the context of weak solutions incorporating the no-slip condition, the following results have 
been obtained. Existence of a weak solution for a FSI problem between a 3 D incompressible, viscous 
fluid and a 2D viscoelastic plate was shown by Chambolle et al. in m, while Grandmont improved 
this result in [23] to hold for a 2D elastic plate. These results were extended to a more general 
geometry in [36] , and to a non-Newtonian shear dependent fluid in [40] ■ In these works existence of 
a weak solution was proved for as long as the elastic boundary does not touch ’’the bottom” (rigid) 
portion of the fluid domain boundary. 

Muha and Canic recently proved the existence of a weak solution to a class of FSI problems 
modeling the flow of an incompressible, viscous, Newtonian fluid flowing through a 2D cylinder 
whose lateral wall was modeled by either the linearly viscoelastic, or by the linearly elastic Koiter 
shell equations [42] . assuming nonlinear coupling at the deformed fluid-structure interface. These 
results were extended by the same authors to a 3D FSI problem involving a cylindrical Koiter shell 
[43] , and to a semi-linear cylindrical Koiter shell [45]. The main novelty in these works was a design 
of a constructive existence proof based on the Lie operator splitting scheme, which has been used 
in numerical simulation of several FSI problems [Ml 13 HU II Hi HD], and has proven to be a robust 
method for a design of constructive existence proofs for an entire class of FSI problems. 

In the present work a further non-trivial extension of the Lie operator splitting scheme is intro¬ 
duced to deal with the Navier slip boundary condition and with the non-zero longitudinal displace¬ 
ment of the structure. Dealing with the Navier slip condition and non-zero longitudinal displacement 
introduces several mathematical difficulties. In contrast with the no-slip boundary condition which 
“transmits” the regularizing mechanism by the viscous fluid dissipation onto the fluid-structure in¬ 
terface, in the Navier slip condition the tangential components of the fluid and structure velocities 
are no longer continuous, and thus information is lost in the tangential direction. As a result, new 
compactness arguments had to be designed in the existence proof to control the tangential velocity 
components at the interface. The compactness arguments are based on Simon’s characterization 
of compactness in L 2 (0,T;R) spaces [53] . and on interpolation of the classical Sobolev spaces with 
real exponents H a (or alternatively Nikolskii spaces 7V S,P [54]). This is new. If we had worked with 
continuous energy estimate, we would have been able to obtain structure regularity in the standard 
space L°°(0, T; H 2 (T)). The time-discretization via operator splitting, however, enabled us to obtain 
an additional estimate in time that is due to the dissipative term in the backward Euler approxima¬ 
tion of the time derivative of the structure velocity. The uniform boundedness of this term enabled 
us to obtain structure regularity in H s (0, T; H 2 (T)), s < 1/2. This was crucial for the existence 
proof. This approach in general brings new information about the time-behavior of weak solutions 
to the elastic structure problems. 

Furthermore, to deal with the non-zero longitudinal displacement and keep the behavior of fluid- 
structure interface “under control”, we had to consider higher-order terms in the structure model 
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given by the bending rigidity of shells. The linearly elastic membrane model was not tractable. Due 
to the non-zero longitudinal displacement additional nonlinearities appear in the problem that track 
the geometric quantities such as the fluid-structure interface surface measure, the interface tangent 
and normal, and the change in the moving fluid domain measure (given by the Jacobian of the ALE 
mapping mapping the moving domain onto a fixed, reference domain). These now appear explicitly 
in the weak formulation of the problem, and cause various difficulties in the analysis. This is one 
of the reasons why our existence result is local in time, i.e., it holds for the time interval (0,T) 
for which we can guarantee that the fluid domain does not degenerate in the sense that the ALE 
mapping remains injective in time as the fluid domain moves, and the Jacobian of the ALE mapping 
remains strictly positive, see Figure [3] Therefore, in this manuscript we prove the existence, locally 
in time, of a weak solution, to a nonlinear moving-boundary problem between an incompressible, 
viscous Newtonian fluid and an elastic shell or plate, satisfying the Navier slip condition at the 
fluid-structure interface, and balance of forces at the fluid-structure interface. 


2 Problem description 

We study the flow of an incompressible, viscous fluid through a 2D fluid domain whose boundary 
contains an elastic, thin structure. The fluid and structure are fully coupled through two coupling 
conditions: the Navier slip boundary condition, and the dynamic coupling condition describing the 
balance of forces at the elastic structure interface. The flow is driven by the data, which includes 
the case of the time-dependent dynamic pressure data prescribed at the “inlet and outlet” portion 
of the fixed boundary, denoted in Figure [2] as Ti and r 3 . The reference fluid domain, denoted by 
C R 2 , is considered to be a polygon with angles less than of equal to 7r, see Figure[ 2 ] If we denote 
by Ti, i = 0,..., to the faces of ft, without the loss of generality we can assume that To is compliant, 
and Ti,..., to is the rigid portion. We denote the rigid portion by E = UI”=i b, and the compliant 
portion by T = r 0 = (0, L). 


r 0 



Figure 2: An example of fluid and structure domains. 

As the fluid flows through the compliant domain, the elastic part of the boundary deforms, giving 
rise to a time-dependent fluid domain which is not known a priori. We denote by tp(t,.) : Cl R 2 , 
t G [0, T) the time-dependent deformation of the fluid domain determined by the interaction between 
fluid flow and the elastic part of the fluid domain boundary. We will be assuming that <p is such that 
the rigid portion E of the fluid domain remains fixed, and that the fluid domain does not degenerate 
in the sense that the elastic structure does not touch any part of the boundary during deformation. 
More precisely, we will be assuming that ip is a C 1 diffeomorphism such that 

y>| S = id and detV<p(i,x) > 0, (i,x) G [0,T) x Cl. 

We denote the displacement of the elastic part of the boundary r 0 by rj{t, x) = <p(f, x) — x, x G r 0 . 
Since To can be identified by the interval (0, L), as mentioned above, 77 is defined as a mapping 
77 : [0,T) x [0,L] R 2 , with 

V(t,z) = {v z (t,z),r) r (t,z)), z G [0, A], 77 ( 0 ) = S z t?(0) = 77 (A) = d z r)(L) = 0. (1) 
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Here r] z and rj r denote the tangential and normal components of displacement with respect to the 
reference configuration To, respectively, and the last set of conditions in (jT|) state that the elastic 
structure is clamped at the points at which it meets the rigid portion of the boundary E. We denote 
by H v (t) = <£>(f, fl) the deformed fluid domain at time t, and by T v {t) = <p(t,r 0 ) the corresponding 
deformed elastic part of the boundary of Vf 1 ( t ). We chose to include r\ as a superscript in the notation 
for the deformed fluid domain and for the deformed elastic structure to emphasize that they both 
depend on one of the unknowns in the problem, which is the structure displacement r\. The following 
notation will be useful in subsequent calculations. The surface element of the deformed structure 
will be denoted by: 

dr 11 = a/(1 + d z r) z (t, z)) 2 + {d r rj r (t, z)) 2 dz = S v (t,z)dz , 

the tangent vector to the deformed structure will be denoted by r v (t, z) = d z ip v (t, z), and the outer 
unit normal on T v (t) at point < p(t, z) will be denoted by i/^t, z). 

The fluid. The fluid flow is governed by the Navier-Stokes equations for an incompressible, 
viscous fluid defined on the family of time-dependent domains Ct^ft): 

M9.U + U VU) = V-„, j infl , (()i(e(0 _ r)i (2) 

where pp denotes the fluid density, u is the fluid velocity, <x = —pi + 2/zD(u) is the fluid Cauchy 
stress tensor, p is the fluid pressure, p is the kinematic viscosity coefficient, and D(u) = ^(Vu+V T u) 
is the symmetrized gradient of u. 

The approach presented in this manuscript can handle different types of boundary conditions 
prescribed on the rigid boundary E. More precisely, on each face Tj of the rigid boundary E we 
prescribe one of the following four types of boundary conditions: 

1. Dynamic pressure: p + |u| 2 = P i7 u • r = 0, prescribed on T, : , i g /, 

2. Velocity (no-slip): u = 0, prescribed on T i: i £ II, 

3. The Navier slip boundary condition: u • u = 0, u • r + c^err/ • r = 0, prescribed on T*, i g III , 

4. The symmetry boundary condition: u • v = 0, d v u T = 0, prescribed on Tj, i g IV, 

where /, II , III and IV denote the subsets of the set of indices {1,..., m} such that the boundary 
condition of type 1, 2, 3 or 4 is satisfied. We note that non-homogeneous boundary conditions can 
also be handled with additional care. 

The problem is supplemented with the initial condition: 

u(0,.) = u o . (3) 


To close the problem, it remains to specify the boundary conditions on the elastic part of the 
boundary T I? (t). For this purpose we introduce the elastodynamics equations modeling the motion 
of the elastic structure, and the two-way coupling between the structure and the fluid motion. 

The structure. The elasto-dynamics of thin structure T I) (t) will be given in terms of dis¬ 
placement with respect to the reference configuration T = To = (0, L) (Lagrangian formulation). 
To include different shell models, we formulate the elastodynamics problem in terms of a general 
continuous, self-adjoint, coercive, linear operator C e , defined on Uq(0, L), for which there exists a 
constant c > 0 such that 

(£ev,v) > c M\ 2 Hi(ry Vj 7ei?o(r)> (4) 

where (■, ■) is the duality pairing between H§ and H~ 2 . The structure elastodynamics problem is 
then given by: 


pshdtt.i 7 = ~£eV + f, 

V{t, 0) = d z r)(t, 0) = r](t, L) = d z r)(t, L ) = 0, 
*7(0, z) =*7 0 , d t r){0,z) = v 0) ; 


*er = (o,L), t g (o,t), 

t g (0,T), 
zer=(0,L). 


(5) 

( 6 ) 
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where ps is the structure density, h the elastic shell thickness, f is linear force density acting on the 
shell, rj = (j)z-,yr) is the shell displacement, and r) 0 and Vo are the initial structure displacement 
and the initial structure velocity, respectively. 

The fluid and structure equations are coupled via the following two sets of coupling conditions: 

• The kinematic coupling condition (Navier slip condition): 

Continuity of normal velocity on T n (t) : 

dtV(t, z ) ■ v v (t, z) = u(ip{t, z)) • v v (t, z ), (t, z) G (0, T ) x T, 

(7) 

The slip condition between the fluid and thin structure on T^ (t) : 
z) - u z))) ■ T v (t, z) 

= aa(4>{t , z))v v (t, z) ■ r\t, z), (t, z ) G (0, T) x T. 

• The dynamic coupling condition: 

Pshd tt r)(t, z) = —C e r\{t, z) - S' n (t,z)cr(<fi(t,z))i' ri (t,z% ( t,z ) G (0,T) x T, (8) 

stating that the structure interface elastodynamics is driven by the jump in the normal stress 
across the interface, where we have assumed, without the loss of generality, that the normal 
stress on the outside of the structure is equal to zero. The term 5^, which multiplies the normal 
fluid stress cru 71 , is the Jacobian of the transformation between the Eulerian and Lagrangian 
formulations of the fluid and structure problems, respectively. 

Notice that there is no pressure contribution in the slip condition (J7| . The pressure contributes only 
through the dynamic coupling condition ([8]). 

In summary, we study the following problem. 


Find (u, rf] such that the following holds: 

The fluid equations: 


Mft u + uVu) = V.„, j in mf)t te(0 , r) . (9) 

The elastic structure (boundary conditions on (0,T) x T): 

Pshd tt r](t, z) + C e rj{t, z) = -S ri (t,z)(r(<p(t,z))v ri (t,z), (10) 

dtr](t,z) ■ v v {t,z) = u{ip{t,z)) ■ ^(M), (11) 

(d t r}{t,z) - u(ip{t,z)) ■ r' ) (t,z) = acr(^{t, z))i/ v (t, z) ■ T v (t, z), (12) 


with 

r?(t, 0) = d z r](t, 0) = rj{t, L) = d z r](, 

Boundary conditions on E: 

p+ ^|u| 2 = Pi, U - T = 0 

u = 0 

u • v = 0, u • t + aicru • r = 0 
u • n = 0, d v u T = 0 

with u T denoting the tangental component of velocity u. 

Initial conditions: 

u(0,.) = u 0 , 77(0,.) = T] 0 , d t 


t,L) — 0, t G (0, T); 


on Tj, sGl, 
on Tj, i G II, 
on Tj, i G III, 
on Tj, i G IV; 


(13) 


-77(0,.) = v 0 . 


(14) 


The initial data must satisfy the following compatibility conditions: 
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• The initial fluid velocity must satisfy: 

u 0 G L 2 { fl 0 ) 2 , V • u 0 = 0, in fl°, 

u 0 • v = 0, on Tj. i £ II U III U IV , (15) 

u 0 • v 0 = v 0 - i^o, onT°, 

where = U^(0), T° = F'(O), u 0 = i^(0, .). 

• The initial domain must be such that there exists a diffeomorphism ip° £ C ,1 (f2) such that 

= fi°, det Vy>° > 0, (y> - I)|r = »7o> ( 16 ) 

and the initial displacement rj 0 is such that 

llr/ollfl-ii/e < c, where c is small. (17) 

We aim at proving the existence of a weak solution to this nonlinear moving boundary problem. 

Before we continue, we note that condition © on the smallness of the H 11 / 6 norm of ri 0 is 
somewhat artificial, and is stated for technical purposes. This condition simplifies the analysis 
presented in this paper which uses Grisvards’s regularity results for elliptic problems on polygonal 
domains fl, which includes our reference domain. With some additional technicalities, we could have 
obtained the same existence result by considering the reference domain to be the initial configuration 
of the fluid domain, which may not be a polygon. In that case we would not need condition © 
but the existence proof would become more technical. 


3 Weak formulation 


3.1 A Formal Energy Inequality 

To motivate the solution spaces for the weak solution of problem ©-© we present here a prelim¬ 
inary version of the formal energy estimate, which shows that for any smooth solution of problem 
|9|-(14) the total energy of the problem is bounded by the data of the problem. The formal energy 
estimate is derived in a standard way, by multiplying equation s pj ) and J 5 | by a solution u and rj. 
respectively, and integrating by parts. The coupling conditions (10) and (121 are used at the bound¬ 


ary T 17 (f), and boundary conditions (13) are used on the fixed portion of the boundary E = UT^ 
We will use u T and ?y T to denote the tangential components of the trace of fluid velocity and of 
displacement at the moving boundary, respectively. We obtain that any smooth solution (u, 77) of 
problem (|9|)-(|l4|) satisfies the following energy estimate: 


1 d 


^ ^(Pi ? ll u lli 2 (n’j(t)) + Psh\\d t ri\\%2(r) 


\'h 2 ( r)) 


|D(u) 


lL 2 (f2>7(i)) 


■f \\ u t - dtV T \\L2(r*(t)) + © — \\ u t\\ 2 l2(y.) < C, 


(18) 


ie/// 


where C depends on the initial and boundary data, and constant c in front of the H 2 - norm of r) is 
associated with the coercivity of the structure operator £ e , see equation Q. 

Before we can define weak solutions to problem © -fll4| we notice that one of the main difficulties 
associated with studying problem © ~( |14[ ) is the moving fluid domain, which is not known a priori. 
To deal with this difficulty a couple of approaches have been proposed in the literature. One 
approach is to reformulate the problem in Lagrangian coordinates (see e.g. [13 Ell). Unfortunately, 
since problem © -(fl4|) is given on a fixed, control volume, with the “inlet” and “outlet” boundary 
data, Lagrangian coordinates cannot be used. In this manuscript we adopt the second classical 
approach and use the so-called Arbitrary Lagrangian Eulerian (ALE) mapping (see e.g. 7, H ii loTTl) 
to transform problem © -(fl4| to the fixed reference domain H. This will introduce additional 
nonlinearities in the problem, which will depend on the ALE mapping. In the next section we 
construct the appropriate ALE mapping and study its regularity properties, which we will need in 
the proof of the main theorem. 
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3.2 Construction and regularity of the ALE mappping 


Definition of the ALE mapping. Motivated by the energy inequality (181 we assume that 
displacement 17 satisfies 


V 


G L°°(0, T; Hq(T) 2 n W 1,oo (0, T; L 2 (r)) 2 C ,o ’ 1 - /3 (0,T; F/f (T)) 2 , 


(19) 


for 0 < /? < 1. The inclusion above is a direct consequence of the standard Hilbert interpolation 
inequalities (see e.g. [31]). We denote the corresponding deformation of the elastic boundary by <p v , 

i.e. 

¥> v (t , z) = id + r](t, z), (t, z) G [0, T] x T. 

We consider a family of ALE mappings A,, parameterized by 77 , 

A„(t) : H —► Q 77 (t), 

defined for each 17 as a harmonic extension of deformation i.e. A r) (t) is defined as the solution 
of the following boundary value problem defined on the reference domain Q: 

AA,(f, •) = 0 in fi, 

A»;(t)| r = <A'(*,•), (20) 

A,,(t)| S = id. 

Regularity of the ALE mapping. Since domain O is a polyhedral domain with maximal 
angle 7r, we can apply Theorem 5.1.3.1 from Grisvard [23] , p. 261, to obtain the following regularity 
of the ALE mapping: 

l|Ar,(i)||w 2 . 3 (n) < C||r 7 (t, .)|| W 5/3,3( r ). 

We can further estimate the right hand-side by the Sobolev Embedding Theorem 


hit, Ollw 6 / 3 ' 3 ^) < C'||77|| ff »( r ), s > 11/6, (21) 

and so 

l|A f; (t)||w 2 '3(n) < C\\ri(t, .)||w s / 3 . 3 (r) < C\\ri(t, -)||jr«(r), s > 11/6. (22) 

By using the Sobolev Embedding Theorem to estimate the W 2,3 -norm of A v (t) we obtain that A v 
has a Holder continuous derivative, namely 


|A„(t)|| CU1/3(0) < C'||r 7 (t, .)||fl-« ( r), a > 11/6. 


(23) 


Finally, we notice that time t is only a parameter in the linear problem (201, thus the regularity 
properties of A, ; with respect to time are the same as the regularity properties of 17 . Now, since 
17 G C'°’ 1 / 12 (0, T; L)) 2 , as shown in (19) for /3 = 11/6, and from the inequality (23) with 

s = 11/6, we have 


1. A, G C°’ 1/12 (0,T;C 1,1/3 (ft)) 2 , and 


2 . 


A f) || c .o,i/i2( OT;C .i,i/3(n)) < C||j7|| C o, i/i2( 0j T ; ifn/6(r)). 


This regularity result is not optimal, but it is sufficient for the remainder of the proof. 
The ALE velocity. The ALE velocity is defined by 


(24) 

(25) 


= Jt A - 

From the regularity of 17 in (19) we see that £ L 2 (0, T; fF 1//2 (fl)) 2 with 

ll w,? llL'»(o,T;j? 1 / 2 (a)) < C'||9(r7|| i oo( 0 T . L 2( r )). 
Furthermore, the following estimate holds: 

ll w, 1 l// 1 (o,T;ff=+ 1 / 2 (f 2 )) < C'||9 t r7|| i 2( 0 ,T;_H=(r)), s < 1. 


(26) 


(27) 

(28) 
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We shall see later in Proposition [7] that the right hand-side of this inequality is, indeed, bounded. 
More precisely, we will show that rj £ if 1 (0, T; H S (T)), s < 1. 

The Jacobian of the ALE mapping. Let us now consider the Jacobian of the ALE mapping 


J 77 ^,.) =detVA„(i). 


(29) 


From the regularity property (24) of A n we have 

JV e C'°- 1 / 12 (0,T;C°’ 1/3 (n)). 


Now, from the compatibility condition (161 we have that the Jacobian at t = 0 satisfies J(0,x) > 
C > 0, x £ Q. By the continuity of the Jacobian as a function of time, this implies the existence of 
a time interval (0 ,T') such that J v is strictly positive on (0,T'), i.e. we have: 


,F >C> 0, on (0, T') x Q. 


(30) 


Injectivity of A v . In order to be able to use the above-constructed ALE mapping to transform 
the problem onto the fixed reference domain, it remains to show that A v is an injection. A sufficient 
condition for the injectivity of r) is given by the following proposition. 


Proposition 1. Let r/ £ L°°(0, T; H 2 (T)) fl W 1,oo (0, L; L 2 (T)) be such that r/(0, .) = r] 0 , and ry 0 
satisfies conditions © and ©■ Then there exists a T" > 0 such that for every t £ [0, T"] the 
ALE mapping A, ( (f) is an injection. 


Proof. First we notice that because of the linearity of problem (20) and the definition of deformation 
ip v , we can write the ALE mapping in the following form: 


A v (t) = id + B J? (t), 


where B,, is the solution of the following boundary value problem: 


AB,,(f,) = 0 in fl, 

Bjj(t)|r = 
B?j(t)|9n\r = 0- 


Now, we see that the regularity of B^ follows in the same way as the regularity of A, ; with analogous 


estimates in terms of rj in the same norms. Therefore, 
things, that 


B, ; satisfies (25), which implies, among other 


sup |VB^I < c(fi). 
n 


We now use Theorem 5.5-1 from [12] (pp. 222), which we state here for completeness: 

Theorem 1 . (Sufficient conditions for preservation of injectivity and orientation [ 12 ]) 

(A) Let ip = id + if : fl C R” —> R n be a mapping differentiable at a point x £ Q. Then: 


|V-01 < 1 —> detV<p > 0. 


(B) Let Ll be a domain in R n . There exists a constant c( fl) > 0 such that any mapping p = 
id + if £ C' 1 (fl, R n ) satisfying 

sup|V-0| < c(H) (31) 

n 

is injective. 

Statement (B) of the theorem says that every domain has as associated constant c(f2) such that 
whenever (31) holds, the mapping ip is injective. 


We use this theorem, together with the assumptions (16) and ( |17| ), to see that there exists a 
T" > 0 such that the ALE mapping A v (t) is an injection for every t £ [0, T"]. 

□ 
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We now take the minimum between T' and T", and call this new time T again, i.e., 

T = min{T , ,T"}, 


(32) 


where T' is determined from the positivity of the Jacobian J 77 , see (301, and T" is determined from 


the injectivity of A rl , see Proposition [lj This new time determines the existence time interval for the 
weak solution. Our existence result will be local in time in the sense that the maximum T for which 
we can show that a solution exists is determined by the time at which the fluid domain degenerates 
in the sense that either the Jacobian of the ALE mapping becomes zero, or the ALE mapping ceases 
to be injective. Examples showing two types of domain degeneration are shown in Figure [3] The 
degeneration of the fluid domain due to the loss of injectivity of A n shown in Figure [3] left can occur 
because the longitudinal displacement is non-zero. The degeneration of the fluid domain shown in 
Figure [3] right, associated with the loss of injectivity of A v and loss of strict positivity of the Jacobian 
J 77 , can occur even if one assumes that the longitudinal displacement of the structure is zero. 


0 


~r, - 



Figure 3: Two ways a fluid domain can degenerate. Left: loss of injectivity of the ALE mapping 
Ajj. Right: loss of injectivity of the ALE mapping A, ; and loss strict positivity of the Jacobian J 77 . 


3.3 The weak ALE formulation 


As mentioned earlier, we will prove the existence of a weak solution to problem <§ -(fl4|) by mapping 
the problem defined on the moving domain Q rl (t) onto a fixed, reference domain Cl, and study the 
transformed problem on Cl. For this purpose we map the functions defined on the moving domain 
CY n (t) onto the reference domain using the ALE mapping introduced above. We will use super-script 
rj to denote that those functions now depend (implicitly) on r/. More precisely, let f be a (scalar or 
vector) function defined on Cl v (t). Then f 77 , defined on Cl, is given by: 

f v {t,z,r) = f(t, A v (t)(z, r)) = f (t,x,y), 

where (x,y) = A n (t)(z,r) £ f l v {t) denotes the coordinates in Cl v (t). Furthermore, we define the 
transformed gradient and divergence operators by 

Vf^t) := (Vf(t)) 77 = Vf r? (VA r) (t)) _1 , V 17 • f’ 7 = tr(V"f’ 7 ). (33) 

It will be useful in the remainder of the paper to obtain a relationship between the gradient V 77 
and symmetrized gradient D 77 ^’^ of the fluid velocity u defined on Cl. This is typically given via 
Korn’s inequality. However, the problem is that the fluid velocity u 77 (t,^,r) for t £ (0,T), which 
is defined on the fixed domain Cl, is coming from a family of velocities u defined on domains II 77 (t) 
for t £ (0, T), via a family of ALE mappings, all depending on r/. We now show that if r) is “nice 
enough”, there exists a uniform Korn’s constant, independent of the family of domains, such that a 
version of Korn’s inequality holds. More precisely, the following Lemma holds true. 

Lemma 1. (The “transformed” Korn’s inequality) Let 

1. rj £ L°°(0, T; H 2 (T)) n W 1 ’ oo (0, L;L 2 (T)), and 

2 . 77(0, .) = rj Q where i) G satisfies conditions © and ©■ 

Then there exists a time T' > 0 and constants C,c > 0 depending only on ||?7||£«>(jj2)nw 1 > 00 (.L 2 ).! 


such that for every u 77 £ H 1 (Cl ) 2 satisfying boundary condition (13) on E, the following transformed 
version of Korn’s inequality holds: 

c||D 77 (‘-)(u)|| i 2 (n) < ||V 77 ^’’- > (u)|| i 2 (n) < C||D 77 ^’'^( U )|| i 2 ( Q ) , t £ [0,T 7 ]. 

The time T' > 0 is determined by the injectivity of the ALE mapping A v . 
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Proof. From Proposition [I] we deduce the existence of a T' > 0 such that A r/ (f) is injective for every 
t £ [0,T 7 ]. Now, the statement of the Lemma follows from the results in [SB] (Lemma 1 and Remark 
6 ) in the same way as in [42]. Namely, for each fixed t £ (0, T') we map u back to the physical domain 
Q r/ (f) and apply Korn’s inequality there in a standard way, using the classical Korn’s constant which 
depends on domain Q v (t) defined by 77. Due to the regularity of 77 given by conditions 1. and 2. 
in the statement of the Lemma, and due to the uniform (in t ) estimate (25), it follows that the set 
{A, ; (f) : t £ (0,T')| is compact in W 1 , 00 (n) 2 , from which the existence of universal Korn constants 
c and C follows. □ 

To define the ALE formulation of problem ([9])-( 14) we recall the definition of the ALE velocity 
given in (26) and define the ALE derivative as a time derivative evaluated on the fixed reference 
domain: 

a t f,o =5 t f+(w". V)f. (34) 

Using the ALE mapping we can rewrite the Navier-Stokes equations in the ALE formulation as 
follows: 

<9t u |n + (u — w 77 ) • Vu = V • er in (35) 

Here, the terms 9 t un and w 77 , which are originally defined on fl, are composed with the inverse of 
the ALE mapping, which maps them back to the moving domain Q v (t). 

Our goal is to define weak solutions to problem |9|-(14) on the fixed, reference domain S2. The 
first step is to introduce the necessary function spaces on fl. For this purpose we notice that the 
incompressibility condition in moving domains V • u = 0 transforms into the following condition on 
D: 

V 17 ■ xT 7 = 0, 

which we will call the transformed divergence-free condition. 

Motivated by the energy inequality (18) we can now define the function spaces associated with 
weak solutions of problem f|)-([T4]). We define the analogue of the classical function space for the 
fluid velocity with the transformed divergence-free condition: 

V v F = {xT' = «,<) € H\n) 2 : V* 7 • u 77 = 0, u 77 • r 77 = 0, on r i; i £ I, 


u 


= 0 on Tj, i £ II, u 77 • u v = 0, on Fj, i £ III U IV}. 


The corresponding space involving time is given by: 


W V F ( 0 , T) = L°°( 0 , T; L 2 (D)) n L 2 ( 0 , T; V v F ). 


The structure function spaces are classical: 


W s ( 0 , T) = TU 1 ’ oo ( 0 , T; L 2 ( 0 , L)) 2 n L°°( 0 , T; H 2 (T)) 2 . 


(36) 


(37) 


The solution space for our problem with the slip boundary condition must incorporate the 
continuity of normal velocities: 


W(0,T) = {(u 77 ,^) € W£(0 ,T) x W S (0,T) 
The corresponding test space is defined by 

Qi(0,T) = {(q 77 , if) £ Cl([0,T);V F x H 2 (T) 2 ) : 


u 


T 


t|r 


v 


= d t r) ■ v n }. 


jV = 


= z) ■ V 11 }. 


(38) 


(39) 


To obtain the weak formulation on U we first consider our problem defined on moving domains 
Q v (t). Take a test function q defined on fU(f) for some 77 , such that the corresponding test function 
q 77 defined on the fixed domain belongs to the test space Q v (0,T), namely, (q' 7 ,i/>) £ Q n (0,T). 
Multiply (35) by q, integrate over U 77 (t), and formally integrate by parts. We obtain the following. 
For the convective term we have: 




((u — w 77 ) • V)u • q = 


1 




((u — w 77 ) • V)u • q — 


1 


/0’7(t) 


((u — w 77 ) • V)q • u 




/ri(i) 


(u - W 77 ) • ^(u ■ q)^ 77 + 


/n»(t) 


(V • w ?7 )u ■ q — 


2 ^ 

iei' 


|uj 2 q v . 
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For the diffusive part of the Navier-Stokes equations we have: 


/ (V ■ <t) ■ q = 2/i / D(u) : D(q) - / 

/Oi(i) dO’ i(t) dooi(t) 


• q, 


where the second term on the right hand side can be expressed as follows: 

q = f f(eriA • v>r] ) c L • 1/1 + ( cri/ n ■ t v ) q • T^dV^ 

drift) ' ' 


crv 1 • q = 


/aoift) 


i&I ' 


PQu 


E /< 


cri/ • rjq • r 


ieiii'- 


I ({(TV 11 ■ V *)ip • I/" + -(<9 t 77 - u) • r^q • t ,? )W’ 7 
drift) ' a ' 


+ £/.»*-£ L v 


U t Q t • 


iei ' 


iein 


/r* 


We sum all the integrals, and take into account the coupling conditions (10)-(12) holding along 


r r/ (t), and boundary conditions ( f3) holding along E. To deal with the non-zero data on T^, i £ I we 
introduce the “source term functional” R which will collect the two terms corresponding to dynamic 
pressure data defined on T,,* £ I: 


(R,q):=J2 Pi Vqe C£([0,T) X fi), 


(40) 


i£l ' 


where v is the unit outward normal to the boundary T;, i £ I. Then we integrate the entire expression 
with respect to time over (0, T) to obtain the following weak ALE formulation of problem 
defined on Q v (t): 

PF L In ( ( StU |n ‘ q+ ' V ) u ' q_ ^ u_ ' V ) q ' u ) + ' w?, ) u ' q ) 

+ / I 2 fj, [ D(u) : D(q) + 'S’' I —u T q T + — I {u T r, - d t g T v)q T vdr ri {t) \ dt ( 41 ) 
do [ dni(t) iein-' 1 '* ai a drift) J 

+Psh [ [ d^rjipdzdt+ [ (C e rj,ip)dt + — / [ {d t g T v — u T n)(p T n S' 1 dzdt = [ (R, q), 
do dr do a do dr do 

for all (q, ip) such that (q’ 7 , ip) £ 2^(0, T), where g T •, = q r v and g vV = qv v denotes the tangential 


and normal component of q, respectively, and the source term functional R is defined in (40) to 
account for the non-zero dynamic pressure boundary condition on T,;, i £ I. 


We now transform (41) to the fixed reference domain via the ALE mapping A r; . To do that we 
first compute the integral involving the ALE derivative d t uj^: 


fj 

I o dni(t) 


5 t u [n • q = 


J^u 17 • q" 


'o d n 
■T 


f 1 [ d t .rn r < • q r) - [ T [ f JoUo-q^O,.). 

do do do do do 


Now, since 

d t r = J v V v • w" (42) 

(see e.g. j|25j, p. 77), the above expression reads 

[ T [ • q = - [ T [ (V 17 • w") (iT 7 • q 17 ) - r u 17 • cW'} - [ Jou 0 • q"(0, .)• 

do doift) do do do 
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Now, we can define the weak solution on the fixed reference domain. 


NOTATION. To simplify notation, from this point on we will omit the superscript rj in iT and 
since everything will be happening only on the fixed, reference domain, and there will be no place 
for confusion. 


Definition 1 . (Weak solution) We say that (u, 77) £ W’ r, (0,T) is a weak solution to problem 
([9]) -(141 defined on the reference domain fl, if for every (q, -0) £ Q^OjT) the following equality 
holds: 


1 0 JTi a i 


u r q T 


(43) 


Y J Q f ^(((u- w?? )-V 17 )u • q - ((u - w”) • V")q • u - (V 17 • w^q • u) 

' J'u-dtq+ [ 2p [ T'D^u) : D ? '(q) + V 
/0 Jn Jo Jn ieII 

+ — / f (u T r, — d t rj T ’i)q T vS ri dzdt - psh [ [ d t rid t xfdzdt + [ (Ceiq,^) 

a Jo Jr Jo Jr Jo 

+ - f f (d t r] T r, - u T v)ipT-nS' 1 dzdt = [ (R, q) + f J 0 u 0 ■ q(0) + [v 0 -ip, 
a Jo Jr Jo Jn Jr 

where J v is the determinant of the Jacobian of ALE mapping, defined in (29), and w ,; is the ALE 
velocity, defined in (26). 

We note that Definition |T] makes sense even in the case when the ALE mapping A v is not 
injective. However, in that case the weak formulation (43) is not equivalent to problem 
even for smooth solutions. 


4 Main Result 


We are now in a position to state the main result of this work. For this purpose we shall assume 
that all the moving domains Ll v (t) are contained in a larger domain f2 max . Indeed, it will be shown 
later, see Corollary [2j that this will be the case. We note that the only reason for this assumption 
is to assume a certain regularity of the source term R, associated with the non-zero boundary data 
on E. Namely, we will be assuming that R £ L 2 (0, T; R 1 (D max )'), where (H 1 )' denotes the dual 
space of H 1 . In the case when only the dynamic pressure boundary data is different from zero, the 
introduction of the source term R is not necessary. However, we state the Main Result in general 
terms, and consider H max to be the union of all the moving domains f V^ft). 


Theorem 2 . (Main result) Let all the parameters in the problem be positive (this includes the fluid 
and structure densities pf and ps, structure thickness h, and the slip-condition friction constants a 
on the moving boundary T v (t), and at’s on the rigid boundary Tjjj). Moreover, let the source term 
functional R, defined in (40), be such that R £ L 2 (0, 00 ; -/^(flmax) 7 ). If the initial data Uo £ L 2 (fl°) 
and r ] 0 £ Hq(T) are such that compatibility conditions ( |15|), (16) and ( |17[ ) are satisfied, then there 
exists a T > 0 and a weak solution (u, rj) to problem (|9j)-(14) defined on (0 ,T), such that the 
following energy estimate is satisfied: 


9 (pf ll u lli“(o,T ; L2(nr,( t ) ) + Ps^||9tr7|| 2 OO ( 0 T;i 2( r )) + c||^|lL°°(o,T ; ff 2 (r))) + Pl|D(u)|| 2 2( 0 , T;i 2(Qr,( t) )) 


+ b 


- dtp 


T\\L 2 (0,T;L 2 (r^(t))) 


+ > . — K 


ielll 


OLt 


ll/ 2 (0,T;L 2 (r\)) < Eo + C 1 |R| 


L 2 (o,T-m(nv(t))') 


(44) 

where C depends only on the initial data and on the parameters in the problem, and Eq is the kinetic 
and elastic energy of the initial data. 
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The proof of this theorem is based on the following approach. We design a partitioned time- 
marching scheme by using the time-discretization via Lie operator spitting. We first separate the 
fluid and structure subproblems and then semi-discretize the resulting sub-problems with respect to 
time. The time interval (0, T ) is subdivided into N sub-intervals (0,1 1 ), • • • , (tjv_i, f jv) each of width 
At, and the fluid and structure sub-problems are solved on each sub-interval. First, the structure 
sub-problem is solved on (L_i, U) using for the initial data the solution of the fluid sub-problem from 
the previous time step, and then the fluid sub-problem is solved on (£j_j, t,) using for the initial 
data the solution of the just calculated structure sub-problem. The process is repeated for each 
sub-interval of (0,T). This defines an approximation of the solution of the coupled FSI problem on 
(0,T). On each sub-interval the transfer of information between the two sub-problems is 

achieved via the “initial data” at L_i. At each time step only one iteration for the fluid sub-problem 
and one for the structure sub-problem are sufficient to obtain a stable and convergent algorithm. The 
goal is to show that as the time-discretization step At —> 0, the sequence of approximate solutions, 
described above, converges to a weak solution of the coupled FSI problem. 

A crucial step in this approach is the way how the fluid and structure sub-problems are designed. 
In particular, we had to be careful to take into account the well-known problems related to the 
so called “added mass effect” [5] to design the fluid sub-problem in such a way that the fluid 
and structure inertia are kept close together via a Robin-type boundary condition for the fluid sub¬ 
problem. In contrast with the no-slip boundary condition, the slip boundary condition is more tricky 
to deal with because of the lack of continuity in the tangential component of the velocity at the 
fluid-structure interface, and so the smoothing of the interface due to the viscous fluid dissipation is 
no longer transferred to the structure in the tangential direction. New compactness arguments based 
on the theorem of Simon [53J and on interpolation of classical Sobolev spaces with real exponents H s 
(or alternatively Nikolskii spaces) will be used to obtain the existence result. Details are presented 
next. 


5 Approximate solutions 

We construct approximate solutions to problem by using the time-discretization via Lie 

operator splitting. 


5.1 Operator splitting scheme 

Let At = T/N be the time-discretization parameter so that the time interval (0, T) is sub-divided 
into N sub-intervals of width At. On each sub-interval we split the problem into a fluid and 
structure sub-problem, and linearize each sub-problem appropriately. Each of the sub-problems will 
be discretized in time using the Backward Euler scheme. 

To perform the Lie splitting we must rewrite problem |9])-(14) as a first-order system in time 


— =AX, t G (0,T), 

X| f=0 = X°, (45) 


where A is an operator on a Hilbert space, such that A can be split into a non-trivial decomposition 
A = Ai + A 2 . For this purpose introduce the substitution v = dt 17 , and rewrite the structure 
acceleration in terms of the first-order derivative of structure velocity. The initial approximation 
of the solution will be the initial data in the problem, namely, u° = u 0 , rf = r) Q , and v° = v 0 . 
For every sub-division of (0, T) containing N £ N sub-intervals, we recursively define the vector of 
unknown approximate solutions 



n+i 

U A t 
n+i 
V At 

n+i 

VAt 


,n = 0,1,..., iV — 1, * = 1,2, 


(46) 
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where i = 1,2 denotes the solution of sub-problems defined by A± or A 2 , respectively. The initial 
condition is given by the initial data in the problem. X° = (uo, Vo, r; 0 ) T . 

A crucial ingredient for the existence proof is that the semi-discretization of the split problem 


be performed in such a way that a semi-discrete version of the energy inequality (44) is preserved 


at every time step. This is associated with successfully dealing with the “added mass effect” Pi I! 42]. 
For this purpose we define a semi-discrete version of the total energy and dissipation at time nAt 
as follows: 


E 


i+i/2 = EF I" jnin ,2 + Psh n\,n+i/2\\2 , / ft rr n-\-ij2\ 


ir / T L ^r ) +(£eVlV / ,r 1 ir /2 )),i = OA-: (47) 

D n At = Atp[ J"D(u£ t ) : D(u2J + ^ + — ||«+ - (<4+11^)- 

Jn 2 ivuM* 


(48) 


Throughout the rest of this section, we keep the time step At fixed, and define the semi-discretized 
fluid and structure sub-problems. To simplify notation, we will omit the subscript At and write 
(u n +S,v"+5 


! V 


n+ 5) instead of {u At 2 ,v A p ,r] n At ^). 


THE STRUCTURE SUB-PROBLEM (Differential formulation): 

n+l/2 _ n \ 

- + C e r 7" +1 =0, 

> on T, 


p s h- 


V 


At 

n+1 _ 


(49) 


T) 


At 


= v «+l/2 


r] n+1 { 0) = d zV n+1 (0) = V n+1 (L) = d z -n n+1 (L) = 0 
The weak formulation is given by: find (v” +1 / 2 , r/" +1 ) € R+T) 4 such that 

“ v " +1 / 2 - v" 


Psh 


V 


At 

n+1 _ 


t/> + (£ e + +1 ,</>) = 0, xjj e [Rq(t)]" 


77 


At 


■ P = 


'* +1/2 -¥>, <p e [l 2 (T)] 2 


(50) 


This problem is similar to the structure sub-problem for the fluid-structure interaction problem 
studied in |42j . where the no-slip boundary condition was considered, and only the radial displace¬ 
ment of the thin structure was assumed to be different from zero. Using the same ideas as in [42] 


one can show that the following existence result and energy estimate hold for problem (50): 


Proposition 2. For each fixed At > 0, problem (50) has a unique solution (v n+ 2 , 7 j ra + 2 ) g 


[-ffg(T)] 2 x [Rg(T)] 2 . Moreover, the solution of problem (50) satisfies the following discrete energy 
inequality: 

(51) 


E 


i«+ 5 
At 


\{Psh ||v" 


2 _ y 


c \\V 


n+i 


n 




<El t , 


where the kinetic energy E At is defined in (47), and c is the coercivity constant defined in Q. 

Proof. The proof of this proposition is similar to the proof in Propositions 1 and 2 in [42] . The 
existence of a unique weak solution follows from the Lax-Milgram Lemma as in [42]. 

The energy inequality is obtained by using v n+1 / 2 in place of the structure velocity test function 

_ n+l/2_ n 

in the first term in the first equation in (50), and by using -— At v in place of the structure velocity 


test function if in the second term in the first equation in (50). After a calculation incorporating 
the equality |(a — b)a\ = |(|a| 2 + |a — 6| 2 — |&| 2 ), one gets 

^l|v ” +1/2 ||| 2 + ^!K +1/2 - v "|| 2 2 + i(£ e r 7 "+ 1 / 2 ,r 7 "+ 1 / 2 ) + \{ C e ( r 7 " +1 / 2 77"), (7 f +1 ' 2 - 77")) 

= ^fh n \\h + l(EeV n ,r t n ). 
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We add the term f n J n \u\ t \ 2 on both sides of the above equality, and use the coercivity property 
of C e to obtain the energy inequality (511. □ 


Remark. We would like to draw the attention of the reader to the fact that the estimate of the 

term II 77"+2 — t /"!! 2 2 was made possible by the fact that we are performing the time-discretization 

. 0 

via operator splitting, and study semi-discretized problems in time. As a result, the approximation 
of the structure velocity, which was semi-disretized using the backward Euler method, gives rise to 
the term II 7^™+2 — 77 n || 2 2 , which corresponds to the well-known numerical dissipation term. By the 

n o _ 

iterative application of inequality (51), the sum with respect to n of the differences I r] n + 2 — rj ra 11 2 2 is 

1__i tl n 


uniformly bounded by a constant, which only depends on the initial kinetic energy and the inlet and 
outlet boundary data, as stated in Proposition [4] below. This estimate gives additional information 
about the behavior in time of the structure displacement, which will be crucial for the compactness 
arguments established in Section [7T~2| 


The structure sub-problem updates the position of the elastic boundary rj n+1 , based on which we 
can now calculate the ALE mapping A n+1 as the harmonic extension of rj n+1 , i.e. A n+1 = id+B” +1 
where B n+1 is defined as the solution of the following boundary value problem: 


AB ' l+1 = 0 in n, B ” +1 = r/ n+1 , B^ 1 = 0. 


(52) 


The corresponding discrete version of the ALE velocity and the Jacobian of the ALE mapping are 
defined by: 

A"+i - A r] 


w" +1 = 


At 


J n+1 = det VA” +1 . 


(53) 


THE FLUID SUB-PROBLEM (Differential formulation): 


Pf- 


A t 


+ pf((u™ — w n+1 ) • V 17 T ')ir 


L n+1 

11 

<1 

43 

n+1 

1 


At 


r 7 ? 


t+1 


• a 

/*(n+l)A£ 
/ nAt 


(u" +1 ,p” +1 ) + R" + i 


R 


in fl 


(54) 


,n +1 _ v n+ 1/2 


At 

(u” +1 - v n+1 ) • 
acr* 1 "^ 1 (u n+1 ,p n+1 )zz' 7 " • t- 


= -S" 

= 0 

= (v" +1 - u n+1 ) 


cr' (u n 1 , p 


on T, 


(55) 


This system is supplemented with boundary conditions (13). 


Before we state the corresponding weak formulation, we introduce the following abbreviations to 
simplify notation: 


V" := V" , cr n := a* , v n := i/» , r n := r ^ , < := u" ■ u n , < := u n ■ r". 


(56) 


We now define the weak solution function space for the fluid sub-problem given in terms of the fluid 
velocity u and its trace v on T as: 

V” = {(u,v) G iL 1 (D) 2 x L 2 (T) 2 : V” • u = 0, (u - v) • u n = 0, 

u • r = 0, on Tj, i G /, u = 0 on T,, i G II, u • u = 0, on Tj, i G III U IV}. 
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The weak formulation is defined as follows: find (u” +1 ,v n+1 ) G V" such that 


Pf 


At 


■ q 


PF 

2 


jn+l _ jn 

At 


,n+1 




u — w 


.n+1 


) • V" +1 )u n+1 


• q 


-((u n - w n+1 ) • V n+1 )q • u n+1 ) + 2/i f J" +1 D(u” +1 ) : D(q) + V [ 

Jn iem^ T 


r, a i 


u? +1 q T 


j « +1 - u" +1 )< 7 t ^+i S n+ 1 dz+p s h 


7 n+l _ v n+l/2 


Jr 

= <R n+1 . 


q). 


At 

(q+) G V 


ip + ~ / ( 


(+ +1 - ++ 1 


)+ 




(57) 

This is obtained by considering (43) and by using formula (42) to express J r, (V 17 • w 1 ') = dtJ v . 
Furthermore, in the third and fourth integral we replaced J" with J n+1 for higher accuracy. This, 
however, does not influence the existence proof. Either choice works well. 

Proposition 3. Let At > 0, and J n+1 > c > 0. Then there exists a unique solution (u ra+1 ,v ra+1 ) 
to the fluid sub-problem, (57). Furthermore, the solution satisfies the following semi-discrete energy 
inequality: 


E 


n+1 

At 




|u n+1 -u"||i, (0) 




n+i _ n+1/2 1 |2 

v v lli 2 (r) 


(58) 


+D1V < E At 


+ 1/2 


+ C'At||R” +1 ||( ff i( n )) / , 


where {H 1 ^ 2))' is the dual space of H 1 ^). 


Proof. The existence proof follows from the Lax-Milgram Lemma and the transformed Korn’s in¬ 
equality stated in Lemma [lj see [12]. 

The semi-discrete energy inequality (58) is a consequence of the fact that we have discretized 
our fluid sub-problem given by (57) so that the discrete version of the geometric conservation law is 


satisfied. More precisely, by taking q 


= u n+1 


in the following two terms in the weak formulation (57): 


pf / r 


At 


- + T 


jn+1 _ jr 

A t 


n+1 


we see that with this kind of discretization a semi-discrete version of the geometric conservation law 
is exactly satisfied , i.e. we have: 


Pf 

2 




n+112 


jn\ u n+l _ 


U+) 


PF 

2 


J n |u 


n |2 


(59) 


Thus, the fluid kinetic energy at time (n + l)Af plus the kinetic energy due to the fluid domain 
motion, is exactly equal to the fluid kinetic energy at time nAt. The two terms on the left hand- 
side in (59) appear on the left hand-side in the energy estimate (58), while the term on the right 
hand-side in ([59| appears on the right hand-side in the energy estimate ([58]) . □ 


Proposition 4. (Uniform semi-discrete energy estimates) Let At > 0 and N = T/At > 0. 

Furthermore, let 2 , E^ 1 , and D At be the kinetic energy and dissipation given by ( ]47| ) and ( |48| ) ; 
respectively. There exists a constant C > 0 independent of At, which depends only on the parameters 
in the problem, on the kinetic energy of the initial data Eq, and on the norm of the right-hand side 
|| R|| l 2 (o T-jjqn)') 0 - e ■ on the boundary data), such that the following estimates hold: 

1. E n At ® < C, EI + 1 < C, for all n = 0,..., N - 1, 

N 

2. J2 D At<C, 


71=1 
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JV-1 r 1 x 

3- ( / ^Afl u A| 1 ~ U Atl“ + ll v A| 1 ^ V At 2 llL 2 (r)+ c l! r 7At 1 ^ ? 7Atllff2( r ) + ||v At 2 —V^ t ||^2(r)^ <C. 

n=0 • /n 

In fact, C = £ i o+C , ||R||^ 2 ( 0 w here C is a constant which depends only on the parameters 

in the problem, and R is the term, defined in (40), coming from the dynamic pressure data. 


Proof. The proof of Proposition [4] follows directly from the energy estimates (51) and [58| ), and from 
the Holder inequality applied to the source term R. More precisely, the statements in Proposition [4] 
are obtained after summing the combined energy estimates (51) and (58) over n = 0, • • • , N — 1, 
and after taking into account that 


N-l 


N-l 


At m 


n+ 1 11 2 


U HHCl)) 




/»n+l At 


n —0 


n —0 


Rll 


nAt 


< C!|R'lll 2 ( 0 ) T;(ff 1 (n))')- 


□ 

This is a crucial estimate which will provide uniform boundedness of approximating solutions to 
problem @-([14]), constructed using our semi-discretized scheme based on Lie splitting. However, 
notice that we have so far only defined the approximate values of our solution at discrete points in 
time, given by nAt. We want to define approximate solutions to be defined at all the points in (0, T). 
For this purpose we define approximate solutions to be the functions which are piece-wise constant on 
each sub-interval ((n— l)Af, nAt], n= 1... of (0, T), such that for t £ ((n— 1) At, nAt], n = l...N, 

u A t(t, ■) = u n At , r] At (t, .) = nl t , v At (i, •) = v£ t , v* At {t, .) = v a 7 5 . (60) 

See Figure |4] We define other approximate quantities in analogous way, i.e. 

Aa t(t, •) = A n At , u At (t ,.) = v At , r At (t ,.) = rl t , w a t(t,.) = w l t , S At (t ,.) = Sl t , j a t(t ,.) = Jl t . 

Our goal is to show that there exists a subsequence of the sequence of approximating solutions 



Figure 4: A sketch of u At . 
defined above, which converges to a weak solution of problem 


6 Convergence of approximate solutions 

6.1 Weak and weak* convergence 

We now focus on the sequences of approximate solutions ( V A t) A t: ( u At)At, and (v A t)At, as At —► 0 
(or, equivalently, as N — ► oo). We first show that these sequences are uniformly bounded, inde¬ 
pendently of At, in the appropriate function norms. The main ingredient in showing the uniform 
estimates will be the results of Proposition [4] 

Proposition 5. The sequence (rj A t)At is uniformly bounded in L°°(0,T-, Hq(T)) 2 . Moreover, there 
exists a T small enough such that A At is an injection, and 

Jlt = detV Al t >°, At > 0,n = 1,... AT. 
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Proof. From Proposition [4] we have that E < C, where C is independent of At. This implies 

ll T 7Aill-L'* J (0,T;i/2(r)) < C. 

To show that A.\ t is injective we again use Theorem 5.5.-1 from T2], stated as Theorem [l] in this 
manuscript, in the analogous way as in Proposition [l] First, we fix At and consider the .ftp-norm 
of the difference between the initial data and the approximate solution at time nAt, to notice that 

INa t ~ %INg(r) < NAtllffg(r) + INolNg(r) < 2< N n=l,...,N. 

Furthermore, we calculate 


IN At - *7olU 2 (r) < INa/ - I 7AtlU 2 (r) = At ^2 II v a* 5 IU 2 (r) < CT, 

i—0 2—0 


where we used that rf\ t = r] Q , and Proposition |4j estimate E^ 2 < C, where C is independent of 
At, to bound 

I Nit — r 7olU 2 (r) < CnAt < CT, n = l,...,N. 

Now, we have uniform bounds for INlt — I 7olU 2 (r) and INlt — ^oINgjr)- Therefore, we can use the 
interpolation inequality for Sobolev spaces (see for example [I], Thm. 4.17, p. 79) to get 


VAt ~ T 7ollff ll /®(r) — 2 CT 1 / 12 , n — 1,... ,N. 


Now from (231 and the construction of the ALE mapping we have 

ll^^-At — I|lc°' 1 / 3 (n) < C\\rj n At \\ H ii /e(r) < |NolN 11 /e(r) +2CT 1 / 12 , 


where I is the identity matrix. 

We want to show that the right hand-side is bounded by a constant which is smaller than or 
equal to the constant c(f2) given by (31) in Theorem [I] (B), which guarantees injectivity of the ALE 
mapping. Indeed, from Proposition [4] we see that C above depends on T through the norms of the 
inlet and outlet data in such a way that C is an increasing function of T. Therefore by choosing T 
small, we can make |Na* ~ T 7olli/ 11 / 6 (r) arbitrarily small for n = 1 ,.... ,1V. Furthermore, by using 
the assumption ( fl7| on the smallness of the initial domain displacement, we see that we can choose 
a T" > 0 small enough, independent of At, such that there exists a constant c > 0 giving 


l|VA^ t - 111 ( 70 , 1 / 3 ( 0 ) < C, 

where c is smaller than or equal to the constant c(f!) from Theorem [l] (B), which implies injectivity 
of the ALE mapping VA^ t . 

Finally, from condition (16) requiring that the “initial” Jacobian of ip° is strictly positive, we see 
that there exists a T' > 0 small, independent of At, such that the Jacobian is strictly positive. 
By taking 

T = min{T , ,T"}, (61) 

we obtain the proof of the proposition. □ 

We note that it is this T, given by Proposition [5j that determines the time interval of existence 
of a weak solution to problem ©-([ill)- since T given in ( |61[ ) is independent of At. The time given 
in (32) is a continuous version of the time given by ([61]). 

Next, we show uniform boundedness of the approximating sequences for the fluid and structure 
velocities. 


Proposition 6. The following statements hold: 

1. (vAt)At, ( v At)At are uniformly bounded in L°°(0,T;L 2 (T)), 

( u At)At is uniformly bounded in L°° (0,T; L 2 (fl)) n L 2 (0,T; /J 1 (fl)). 
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Proof. The proof follows from Propositions [4] and [5] and from the uniform Korn’s inequality stated in 
Lemmajl] Namely, statement 2 of Proposition|4]implies that D' ,At (uAt) is bounded in L 2 (0, T; L 2 (VL)). 

Now, we can use the uniform bound for ry At and the uniform Korn’s inequality from Lemma |T] to 
finish the proof. □ 

Lemma 2. (Weak and weak* convergence results) There exist subsequences ( v a*)a t, ( v At)At> 

and (u At ) At , and the functions p G L°°(0, T; Hq(T)), v € L°°(0, T; L 2 (r)), v* G L°°(0, T; L 2 (r)), 
and u G L°°(0, T; L 2 (fl)) n L 2 (0, T; H 1 ^)), such that 

Va t - 1 ‘ V weakly* in L°°(0, T-, Hq(T)), 
v At — *■ v weakly* in L°°(0, T; L 2 (T)), 

v At v* weakly* in L°°(0,T;L 2 (r)), (62) 

u At —^ u weakly* in L°°(0, T; L 2 (Q)), 
u At — ^ u weakly in L 2 (0, T; H 1 ^)). 

Furthermore, 

v = v*. (63) 

Since our problem is nonlinear, we need strong convergence of approximating sub-sequences to 
be able to pass to the limit, and show that the limiting functions satisfy the weak formulation of 
the problem. For this purpose we need a compactness result, which we present next. 


6.2 Compactness 

Compactness arguments in the case when the boundary condition on r(t) is the slip condition is 
different from the compactness argument for the problem in which the boundary condition on T(t) 
is no-slip. In the no-slip case the viscous dissipation of the fluid smooths out the interface providing 
spatial regularity of the interface velocity v At , which is no longer available in the slip condition 
case, since the fluid and structure velocities are no longer equal at the fluid-structure interface. We 
will have some help from the viscous fluid dissipation in the estimates for the normal component of 
structure velocity, but the tangential component no longer “feels” fluid dissipation as before. It is 
because of this that we need to construct a different compactness argument, which we present next. 

To investigate compactness (regularization) in time we introduce the translation in time by h of 
a function /, denoted by Th , as: 

T h f(t, ■) = f(t - h ,.), h G K. (64) 

The following estimates hold for our approximate solutions as they are shifted in time by h: 
Lemma 3. There exists a constant C > 0, independent of At, such that for every h > 0 we have: 

||I)iU At — u A t|| L 2( ftiT;L 2(n)) < C\/h, 

\\T h VAt — v At||L2(/i,T;L 2 ( r )) < Cy/h, 
ll^ v At - v AtlU 2 (h,T;L 2 (r)) < CVh, 

\\ThflAt ~ r 7AtlU 2 (/i,T ; ff 2 (r)) < CVh, At > 0. 

Proof. The proof is analogous to the proof of Theorem 2 in [42]. A summary of the main steps is 
the following. We focus on the first statement given in terms of u At , while the proofs for the other 
two are analoguous. First, from Proposition [4] we immediately have: 

|I^AtU At — UAi || z. 2 ( A *,'T;i 2 (G2)) = Dl'«t 1 -'A.llv<n) a *S CA ‘- 

n 

This implies that the first estimate in the above Lemma holds for “the diagonal” terms for which 
the translation is performed by h that is exactly equal to At. However, we would like to prove the 
statement for an arbitrary translation by h > 0. Let us fix At > 0 and consider the following two 
cases: 0 < h < At and 0 < At < h. We obtain the desired estimates by calculating the following. 
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1. For 0 < h < At we have 


\\ThUAt 


UAi 


L 2 (7i,T;L 2 (fi)) 


E 


pnAt 
I nAt—h 


l|u 


71+1 

At 


x At\ 


2 

L 2 (Q) 


= hJ2\\uTt 1 -u n At\\h { n)<Ch. 

n 

2. For 0 < At < h we write h = l At + s for some l £ N and 0 < s < At, and use the triangle 
inequality to get 


i+i 

\\T h u At - u At \\ 2 L 2 {KT . L 2 m < At ^2 E H u At* _1 - u At*lli 2 r < C(l + l)At < Ch. 

n 7—1 


□ 

In what follows, it will be useful to introduce a slightly different set of approximate functions 
for u, v, and r/ by extending the values of those functions at points nAt to the time sub-interval 
[(n — l)At, nAt] not in a piece-wise constant fashion as before, but linearly. Namely, for each fixed 
At, define u^r, fl At and VAt to be continuous, linear on each sub-interval [( n — l)At, nAt], and such 
that 

u A t(nAt,.) = u A ((nAt,.), v At (nAt ,.) = v At (nAt, .), fi At (nAt,.) = r] At (nAt, .), (65) 

where n = 0,... N. We now observe that 


«ra+l — „n 

r) A t (t) = - -(t — nAt) + rj n , t £ [nAt, (n + l)At), n = 0,..., N — 1, (66) 


dtV At(t) = 


jyll+l _ yjU r] n + 1/2 _ yjTl 


= v n+2 , t e (nAt, (n + l)At), 


At At 

and so, since v At was dehned in (601 as a piece-wise constant function defined via v At (t, •) = v n+ ^, 
for t £ (nAt, (n + 1)At], we see that 


d t f) At = v* At a.e. on (0,T). 


(67) 


The following Lemma will be crucial for establishing compactness of the approximate sequence of 
solutions f] At - 

Lemma 4. There exists a constant C > 0, independent of At, such that for every h > 0 we have: 

\\ T hV A t ~ V A t\\L 2 (h,T-H 2 (r)) < CVh. 

Proof. In the same way as in Lemma [3] we consider two separate cases: 0 < h < At and At < h. 
Case 1: 0 < h < At. We use ( |66| ) and explicitly calculate the straight lines defining the function 
f] At and its translation in time by h to the right, to obtain: 

W^hVAt ~ V A t\\H 2 (r) = J II ^')llff 2 (r) dt 


N-l 


( pnAt+h i 

+ E { / (x72llM»7lt-<t 1 ) + +At-t)(T 7 ^ 1 -2r ? l t + r ? ^ 1 )||^ (r) ) 

n=l dnAt 


dt 


r*(n+l)Ai n +1 _ 

II VAt VAt , 


(— ^)llff 2 (r+} ^ + ff-j) E NaI 1 ~ T 7Atllff 2 (r) ^ Ch. 


N -1 


' nAt-\-h 


At 


n —0 


The last inequality follows from Proposition [4] and from the fact that h < At. 
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Case 2: At < h. We notice that T/f) At = Thrj At and use the following identity (see e.g. [55], p. 
328) 

At N ~ 1 

WVa t ~ ^Atlli 2 (0,T;iT 2 (0,.L)) — ^ II^At 1 “ T ?Atllj/ 2 (0,L)- 

n —0 

From Lemma [3] we have: 

\\ThVAt — VAt \\ L 2 (h,T;H 2 (T)) < W^hVAt — T h T] At ||l 2 (0,T;H 2 (0,L)) 

+ ll^ 1 /i ? 7At — V At II L 2 (0,T;// 2 (0,L)) + \\VAt ~ VAt l|i 2 (0,T;ff 2 (0,L)) < C(At + h + At) < Ch. 

The statement of the Lemma immediately follows by taking the square root on both sides of the 
inequality. □ 

Proposition 7. The following statements hold: 

1. ( u At)At>o uniformly bounded in H s (0,T; L 2 (fl)), 0 < s < 1/2. 

2. (f]At)^t>o Is uniformly bounded in H s (0 7 T; H 2 (T)) fl i7 s+1 (0, T; L 2 {T)), 0 < s < 1/2. 

3. (v At ) A i>o Is uniformly bounded in H s (0,T; L 2 (T)) H L 2 (0, T; H 2 s (T)), 0 < s < 1/2. 

Proof. Since we have already proved that ua* is uniformly bounded in L 2 (0, T; L 2 (fl)), it only 
remains to prove that the following semi-norm is finite: 


UAt|lfl'»( 0 ,T;L 2 (n)) 


|u At (t) - u At (r)|| 2 2(n) 


.0 Jo \t~T\l+ 2 * 

By a simple change of variables h = t — r and Lemma [3] we get: 

f T dh 


dtdr. 


|U A t||#s( 0 ,T;L 2 (Q)) 


' —T N 


l+2s 


u At (t - h) - u At (t)||| 2(f2) dr < 


\h\dh 


— T 1^1 


l+2s ' 


This integral is finite for s < 1/2 and therefore we have proved the first statement. 

The first part of the second statement and the first part of the third statement, i.e. (f) At )At>o and 
(v At )Ai>o are uniformly bounded in H s (0 7 T;H 2 (T)) and H s (0,T;L 2 (T)), 0 < s < 1/2, respectively, 
are proved analogously by using Lemma [3] again. 

Let us now prove the boundedness of (v At ) At >o in L 2 (0,T; H 2 s (T)). First notice that from ( |67| ) 
and the uniform boundedness of (v^ ( ) A i >0 in H s (0 7 T; L 2 (T)), s < 1/2 that we just proved, we see 
that (d t f] At ) is uniformly bounded in H s (0,T- L 2 (T)), s < 1/2. Therefore, 

(t) a t) is uniformly bounded in H s+ 1 (0,T - 7 L 2 (T)) n H s (0, T; H 2 (T)), s < 1/2. 

Now, by the interpolation property (see e.g. [35] Section 1.9.4 p.p. 47) we obtain that ( f ) At ) is 
uniformly bounded in H 1 (0,T-, H 2 s (T)). By using (67) again we conclude that (v^ t ) At>0 = ( d t f)At ) 
is uniformly bounded in L 2 (0, T; iL 2s (r)), s < 1/2. □ 

Notice that the proof of this Proposition heavily relies on the definition of the new approximate 
solution sequences (651 and their properties (66) and (67). In particular, (67) allowed us to obtain 
information about the regularity properties of (v^ t ) At>0 via the regularity of (d t f) At ). It is because 
of this result that we introduced the new definition of approximate solutions given by (65). We 
shall see below that the limits of approximate sequences as At —> 0 do not depend on the type 
of extension of the approximate values of the solution at points nAf onto the time sub-interval 
(nAt, (n + l)At). Therefore, the introduction of the new approximate sequences in (65) to obtain 
additional information about the regularity of the limiting solution is justified. 

We note here a side remark that we could have proved Proposition [7] by using a slightly different 
but equivalent approach, relying on Nikolskii spaces. More precisely, from Lemma[3]we can directly 
conclude that sequences (u )At, (v) At , (v)^ t , (■ r))At. are uniformly bounded in the following Nikolskii 
spaces !V 1 / 2 ’ 2 ((0, T); L 2 (fl)), 7V 1 / 2,2 ((0, T); L 2 (r)) and A 1 / 2,2 ((0, T); H 2 (T)), respectively (see e.g. 
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jo~Tj for the definition of the Nikolskii spaces). Similarly, Lemma [4] gives uniform boundedness of 
(f)) At in A rl / 2,2 ((0, T); H 2 (T)). Proposition [ 7 ] then follows directly from the embeddings of Nikolskii 
spaces into H s spaces [54]. 

We are now ready to state our main compactness result. It relies on the following compactness 
theorem by Simon, stated in 231 as Corollary 5. 

Theorem 3. (Corollary 5 [55] ) Assume that X, B, and Y are Banach spaces, and X C B C Y 
with compact embedding X CC B. Let 1 < p < 00 and 1 < r < 00 . Let F be bounded in 
L p (0,T;X) fl W s ’ r {0,T; Y), where s > 0 if r > p, and where s > 1/r — 1/p if r < p. Then F is 
relatively compact in L p (0,T\ B) (and in C(0,T; B) if p = 00 ). 

Theorem 4. (Compactness) Sets {u At : At > 0} and {v* At : At > 0} are relatively compact in 
L 2 (0,T; H 2 s (Ll)), s < 1/2, and L 2 (0,T-, H 2 s (T)), respectively. 

Proof. The proof of this theorem follows from the estimates obtained in Propositions [6] and [7] and 
by applying Theorem [3] with p = r = 2 and s > 0. □ 

We remark that in contrast with the no-slip condition case studied in |42j . where we obtained 
partial regularity of dtp from the trace of the fluid velocity on the interface and the estimates related 
to the fluid viscous dissipation, here we had to calculate directly all the time-shifts for rj At ,u At , 
and v At , and v At to obtain uniform boundedness in the H s spaces which, combined with the 
Simon’s Corollary 5 and interpolation of classical Sobolev spaces H s with real exponents s, provide 
compactness. 

The compactness result stated in Theorem [4] implies the following strong convergence results. 
Corollary 1. We have the following strong convergence results as At —> 0: 

1. u At -> u in L 2 {Q,T]H 2 s (Ll)), s < 1/2, 

2. v* At -> v in L 2 (0,T-H 2 s (T)), s < 1/2, 

3. v At -> v in L 2 (0, T; H 2 s (T)), s < 1/2. 

To get strong convergence results for the structure displacements 17 At we proceed in the same 
way as in [52|. Namely, from Propositions [ 5 ] and [bj and from di) At = v At , we obtain that r) At is 
uniformly bounded in L°°(0, T; F[q (r)) 2 fl kW’ 00 ^, T; L 2 (P)) 2 . From the continuous embedding 

L°°(0, T ; ILg (r)) 2 ft VF 1 ’ oo (0, T; L 2 (r)) 2 -4 C°’ 1_a ([0, T], H 2 a {T)), 0 < a < 1, 

we obtain uniform boundedness of r) A( in C o,1_Q ([0,T], 14 2Q (r)). Now, to get compactness in 
space we recall that H 2a is continuously embedded into H 2a ~ € . By the Arzela-Ascoli theorem this 
embedding is compact. In fact, by the application of the Arzela-Ascoli theorem to the functions 
in C° :1_ “([0, T], H 2 a (T)) we get the compactness is time as well. More precisely, we obtain the 
existence of a subsequence, which we denote by i) At again, such that 

VAt -> V in C([0, T]; H 2 s (T)), 0 < s < 1. 

Since sequences r) At and 17 At have the same limit i) = 17 G C([0, T]; L[ 2 s (T)), where 17 is the weak* 
limit discussed in Lemma [2] we obtain 

VAt^V in C([0,T];ff 2 s (r)), 0 < s < 1. 

By combining this statement with the continuity in time of 17 , as was done in Lemma 3 of [42] . we 
obtain the following strong convergence results for the structure: 

Theorem 5. The following strong convergence results hold as At —> 0: 

1. r, At ^r,mL°°(0,T;H 2 s (T)), s < 1, 

2. T At 7 lAt 17 in L°°(0,T; H 2 s (T)), s < 1. 
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To pass to the limit in the weak formulation, we need uniform convergence of r] At . This is where 
the fact that we work in 2D rather than in 3D comes into play. Namely, since in our 2D problem T 
is a ID domain, we have that H a (T) is embedded into C 1 (T) for s > 3/2, and so the first statement 
of Theorem [H] implies rj At —► rj in L°°(0, T; C' 1 (T)). Moreover, we have the following Corollary. 

Corollary 2. The following uniform convergence results hold as At —> 0: 

1. r] At —> r] in Z/°°(0, T ; C 1 ^)), 

2. T At n At ^ r, in L^^T-C 1 (7)). 

The second statement in this corollary can be proved using the same arguments as those following 
statement (76) in [32] ■ 

By using this corollary, and the explicit formulas for the normals u At , the tangents t A t , and the 
quantities associated with the ALE mappings A At , one can see that the following strong convergence 
results hold: 

Corollary 3. The following strong convergence results hold for the geometric quantities associated 
with the change of the fluid domain Cl(t): 

1 . i/ At —► */" in L°°(0,T;C(r)), 

2 . T At ^r’MnL oo (0,T;C(r)), 

3 . WAt-^w^in 

I S At ->• S" in L°°(0,T; C(7)), 

5. J A t~^ J v ™L°°(0,T-,C(n)), 

6 . (VAaO" 1 -> (VA,,)" 1 in L°°(0, T; C(f2)). 

We remark that the results of Corollary [3] were not necessary in our previous work [32] because 
only the normal component of structure displacement was assumed to be non-zero. In the current 
manuscript both the normal and tangential structure displacements are considered to be non-zero, 
which introduces additional complications in tracking the change in the measure of the interface 
“surface” deformation that we did not have to deal with before. 


7 The limiting problem 

7.1 Construction of suitable test functions 


Now that we have the strong convergence results above, we are ready to show that the limits, as 
At —> 0, of approximate solutions satisfy the weak form (43) of problem <n> -fll4|). Unfortunately, due 
to the fact that we mapped our problem defined on the moving domain Q(t) onto a fixed, reference 
domain ft, introduces additional difficulties. More precisely, the velocity test functions in the weak 
formulation of the fluid sub-problem (57) now depend of At via their dependence on rj At . This is 
because of the requirement that the transformed divergence-free condition V 17 ** • q = 0 must be 
satisfied. Passing to the limit in the weak formulation of the fluid sub-problem (57) when both the 
test functions and the unknown functions depend on At is tricky, and special care needs to be taken 
to deal with this issue. 

Our strategy is to restrict ourselves to a dense subset, call it X v (0,T), of the space of all 
test functions Q v (0,T), and for every q S X v (0 ,T) construct a sequence of test functions for 
the approximate problems, call them qAt, such that q A t —> q in suitable norms. This approach was 
used in |4'2j for the FSI problem with the no-slip condition and only radial structural displacements, 
see also mmm- Since, here the test space is different because of the slip boundary condition, 
the construction of such test functions is somewhat different. 
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First let us define the domain which contains all the approximate domains 


!)max= |J W*. 

At>0,nSN 


(68) 


Notice that Corollary [2] implies C Umax, t G [0, T], and E C <9f2 max . We define 

dfmax ={r € ^([0,3-); C' 2 (O max )) 2 : V • r = 0,r • r = 0, on T t , i G I, r = 0 on T z , * G //, 
r • v = 0, on r i; i G III U IV}. 

*”(0, T) ={(q, 0) : q (t, •) = r(t, -)|o^(t) ° A„(i), r G A max , (r| r „ - ij>) ■ = 0, 0 G H$ (r)}. 

From the construction it is immediate that X v (0,T) is dense in Q v (0,T). 

We now want to construct the test functions qAt and 0 At for the approximate problems, such 
that qAt -t q and 0 At — > 0 in a suitable space. For this purpose let us fix (q, 0) G X v (0,T), 
q(t, .) = r(t, ,)| Qv(t) ° A , ; (t), r G A max , and define (qAt,0At) to be piece-wise constant in time so 
that: 

qAt(C •) = q r At : = r(nAf, ■)|n’?At(t) ° A n At (t), 


0At(0 = 0At : = 0(«Ai), 


, t G ((n — l)Af,nAf]. 


(69) 


Note that (qAt(t, -),0At(f, •)) e V^t, t G ((n — l)At,nAf], Now, using ideas from [12] we can prove 
the following lemma. 

Lemma 5. For ererj/ (q, 0) G X v (0,T) we have 

(qAt, 0 A t) (q,0) in F°°(0, T ; C 1 (H)) 2 x L°°(0,T; C 1 ^)) 2 . 

We will also need information about the convergence of approximations of d t q, which we define 


by: 


n” +2 - n" +1 

dq A t(t,.):= At At , i€((n-l)Ai,nAi]. 


(70) 


Lemma 6. Fef (q, 0) G (0, T), and let dq At &e defined by (70). Then dq At —> d t qin F 2 (0, T; F 2 (S1)). 
Proof. By the Mean-Value Theorem we have 

= ^(r((n + 2)At,A’A+ 2 (r,z)) - ((n + l)Af, A^f(r, z))) 

= ^(r((n + 2)Ai,Aft 2 (r,*)) ~r((n + l)At,A n A f(r,z)) + 
r((n + 1) At, A£+ 2 (r, z)) - ((n + l)Ai, A^(r, z))) 

= $r((n + 1 + /3)At,A”+ 2 (r,z))) + Vr((n + l)At, C) At 1 At ( ’ J , 


A"+2 a n+1 
A At A At 


n+2 


where £ = A”+ (r, z) + 7(A^+ (r, z)- A^+ (r, z)),/3, 7 G [0,1]. Notice that-—-- « At , 

and this term is associated with T_AtWAt, which converges strongly to w I) in F 2 (0, T; F 2 (fi)). 
Therefore we have 

dq A t -t <9 t r + Vr • w v = d t q in F 2 (0, T; F 2 (fi)). 


□ 
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7.2 Approximate equations 

We have so far introduced the weak formulation of the coupled problem at the continuous level, and 
have split the coupled problem into the fluid and structure sub-problems. We then semi-discretized 
the two sub-problems, and introduced the semi-discrete weak formulations of those approximate fluid 
and structure sub-problems. What needs to be done next is to define the semi-discrete weak formu¬ 
lation of the coupled problem which approximates the weak formulation of the coupled continuous 
problem. To do that we take the constructed approximate test functions (qAt, '•/’At); multiply them 
by At, and replace the test functions q and ip in the weak formulations for the approximate structure 


and fluid sub-problems (501 and (57) with the test functions At(q At , ip At ). We add the two weak 
formulations together, and sum w.r.t. n = 0,..., N — 1. The approximating solutions (rL At ,r/ At ) 
satisfy the following variational form of the semi-discretized (approximate) coupled problem: 


Pf 


~ Pf 

TAtJAtdtU At • qAt H-W" 


to Jn 


/ o Jn 


J At — T At JAt 
A T 


-UA t ■ qAt 


+ [ [ J At(((7AtU A t - w At ) • V^^UAt • q At -((T At u At - w At ) • V^^qAt • u At ) 

J 0 J n 

+2 p [ [ JAiD^UAt) : D* 7 **(q A t) + - [ /((u At - v At ) • T At )(qAt ■ r At )S At dz 

Jo Jn a Jo Jr 

+ - / ((VAt - u A t) • F At )(ip At ■ T At )S At dz+ V' / — (u A t • r)(q At • r) 

a J o J r TTlrJr^i 


(71) 


i£lll‘ 


+Psh 


d t v At ip At + {CeV At , ip a t) = (R-At, qAt), (q, ip) e ^(o, r). 


We want to pass to the limit as At —> 0 and show that the limiting functions satisfy the weak 
formulation of problem ([9j)-( 14), given in ( |43[ ). Indeed, by using the convergence results for the 
approximate solutions given by Theorem |4j and Corollaries [2land[3j and by using the convergence 
results for the corresponding test functions, given by Lemma[5j we can pass to the limit directly in 
all the terms, expect the ones associated with the geometric conservation law of the ALE mapping, 


i.e. the first two terms in the first line in (71). 


7.3 Discrete v.s. continuous geometric conservation law 

We show that the terms associated with the semi-discrete approximation of the geometric conserva¬ 
tion law associated with our family of ALE mappings A At , converge, as At —>• 0, to the corresponding 
terms associates with the geometric conservation law satisfied by the ALE mapping A 11 appearing 
in the continuous weak formulation ( [43] ). More precisely, we have the following result. 

Proposition 8. For every (q ,ip) € A 7) (0,T) the following convergence result holds 


Pf 


~ Pf 

TAt.JAtdtb A t • qAt H— 


JAt. — TAtJ At _ 

A T 


-UAt • qAt -> 


-pF [ [ J v U ■ dpi - • w v )u • q f j 0 u 0 q(0,.). 

Jo Jn z j q Jn Jq 

Proof. Let us first consider the term that contains the fluid acceleration dtV. A t- We use the definition 
of approximate solutions and test functions, and the summation by parts formula to obtain 



TAtJAtdtU At ■ qAt 


N-l 

E 

n—0 


[ JU u; 

Jn 


n+1 
At 


— U 


At, 


■ = 


[ J N . n JV+l 
/ J At u At UAt 
Jn 
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c 1 r 

- / Jit uL • q L - E / u a7 • - -Ml 1 )- 

•' S i n=0 ^ 


Notice that by construction we have q^ t = 0, n > N. By adding and subtracting q^{ we can 

write 


JV-l Af-l JV-1 . 

E / «S 1 -W 1 q3J- a -J2 t q3J 1 ) = E / <7 W^HE / <Mq£f-qAf K 

n=0 n=0 n=0 


n+1 
At ‘ 


By plugging this calculation back into the above formula we get 


/o Jn 


l\tJ aA^av TAt — — / -^At u At'qAt — 


/O ,/fi 


<7a* ~ TAt JAt _ 

aF 


-UAfqAt- 


?At < 7 At UAt • dqAt • 


/o -/f2 


Therefore, the two terms on the left hand side in the statement of Proposition [8] are equal to 


Pf 


= ~Pf 


TAtJAt^i^At ' QAt + 


Pf 


JAt ~ At JAt 


/o Jn 

r-T r 


lo Jn 


At 


-UAt • qAt 


/o Jn 


TAtJAt^-At • *lAt Pf I </o u o • qXt — 


'o Jn 


JAt — TAtJAt 

A T 


-UAt • qAt- 


(72) 


We can pass to the limit in the first two terms on the right by using Lemma [bj Corollary [3j and 
Corollary [l] Passing to the limit in the third term on the right is not that straight forward. We 
want to show that 


/7 


JAt ~ TAtJ At _ 

At 


-u At ■ q -t 


fj 


J 17 (V^ 


>-q. 


The main problem is passing to the limit in the term (JAt — TAtJ At) I At, since it is not clear that 
it converges to J^, which, as stated in (42), is equal to d t J v = J V (S7 V ■ w v ). 

The plan is to explicitly calculate Jai ~ TAtJ At by using the mean value theorem. For this 
purpose we first recall that the integral with respect to time of Ja* — TAt JAt is equal to the sum 
over n = 1,..., N of the differences J^j) — Jit times At. Now, since 


Jit = detVA^ f > 0, At > 0,n= 


we have 

Jit 1 - JAt = det VA^j 1 - det VA£ t . 

We apply the mean value theorem on the determinant function when the difference VA^j 1 — VA^ t 
is small. The mean value theorem says 

det VA^j 1 - det VA^ ( = D(det)(VA”f )(VA^+ 1 - VA£ t ), (73) 

where D(det)(VA^) denotes the derivative of the determinant function evaluated at an interme¬ 
diate point 


VA”f := V A At + ^(VA”! 1 - VA£ t ), for some /3 £ 


l At ~ v-ri-At)t iur some p t [0,1]. 

The functional D(det)(VA^’f) acts on the difference (VA^ 1 — VA^ t ). 

To explicitly calculate the right hand-side of (73), we use the formula for the derivative of the 
determinant, evaluated at F, acting on U, given by 


D(det)(F)U = (detF)tr(UF- 1 ). 


By using this formula we get 

Jit 1 ~ J At = det (VA"f )tr((VA Af n+1 - VA At n )(VA”f )" 1 ), for some /3 £ [0,1], 
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Now, from (33) we have that the factor on the right hand side containing the trace is equal to 


tr V 


A»+l _ \n \ \ A n +1 _ A n 

^-a t t j (^7A n ’P | ^7 n >@ . a ai _ a a t \/ n P 


At 


At 


n+l 
/ A t > 


where V n/3 = V 17 " + /? (v 77 " 
we get 


- V 17 ' 


and by denoting detVA^’f with in the spirit of (531, 


rn-t-l __ t 
J At J 


At 


At 


= J2f (v^-w"! 1 ) 


Thus, we have just calculated that 


jn+1 _ jn 

J At J At „n+1 

At Ua ‘ 


( ilV = / r 




w Ii+1 lii n+1 
W A t ) U At 


QaI 1 - 


By taking the sum over n = 0,..., N — 1 in equation (|74|) we get 
JAt — TAt JAt 


rT-At 


pT — At p 

/ / • / At( V ’ ?A ‘ ' w At)u A f ' qAt- 

Jo J a 


At 


-UAt ■ qAt = 


Now we can pass to the limit as At —> 0 to obtain: 

f T f JAt ~ TAt JAt 


At 


-u At • q -t 




>-q. 


(74) 


(75) 


(76) 


Finally, with this conclusion we can pass to the limit in (72) to obtain 


Pf / TAtJAtdt^iAt • qAt + ~ 


JAt — TAt JAt _ 

aT 


-UAt • qAt 


-Pf [ [ J v u d t q- [ J 0 u o q(0, .)-^f [ [ J n (V n ■ w r ')u 
Jo Jn Jn 1 Jo Jn 


■ q, 


which is exactly the statement of the Proposition. □ 

Therefore, we have shown that in the limit as At —> 0, the approximate solutions constructed 
in Section [5] based on the Lie operator splitting scheme, converge to a weak solution of problem 
0 -(fl4|. More precisely, we have shown the following result. 

Lemma 7. There exists a T > 0, and a subsequence of approximate solutions (uAt, VAt)> constructed 
in Section^ such that (uAt,»7At) converges, as At —>■ 0 , to a function (u, if) £ W v (0,T), which 
is a weak solution to problem § -( |14| ) in the sense of Definition^ The weak form in Definition^ 7] 
holds for all the test functions (q, if) £ X v (0 ,T) which are dense in Q v (fi,T), and are obtained as 
the limits of the test functions (qAtj^At) £ Vaj constructed in SectionVLl I 


We are now ready to complete the proof of the main existence result, stated in Theorem [2] From 
Lemma[7]we obtain the existence of a weak solution defined on the time interval (0, T), where T > 0 
is determined by (61). To obtain the energy estimate ( [44] ) from Theorem [2] we consider discrete 
energy inequalities stated in points 1. and 2. in Proposition [4] and let At —> 0. Due to the lower 
semi-continuity property of norms, we can take the limit in points 1. and 2. in Proposition [4] to 
recover the energy estimate (441. 

This concludes the constructive proof to the main existence result, stated in Theorem [2] 


8 Conclusions 

This paper provides a constructive existence proof for a weak solution to a nonlinear moving bound¬ 
ary problem between an incompressible, viscous fluid and an elastic shell, with the Navier slip 
boundary condition holding at the fluid-structure interface. Due to different types of boundary 
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conditions holding at each piece of the fluid domain boundary the usual vorticity formulation, com¬ 
monly used in a good agreement with the Navier slip boundary condition, does not appear helpful 
for this problem. The problem is motivated by studying fluid-structure interaction between blood 
flow and cardiovascular tissue, whether natural or bio-artificial, which include cell-seeded tissue con¬ 
structs, which consists of grooves in tissue scaffolds that are lined with cells giving rise to “rough” 
fluid-structure interfaces. To filter out the small scales of the rough fluid domain boundary, effective 
boundary conditions based on the Navier slip condition have been used in various applications, see 
a review paper by Mikelic m and the references therein. The present work is the first existence 
result involving the Navier slip boundary condition for a fluid-structure interaction problem with 
elastic structures. Dealing with the slip condition introduces several mathematical difficulties. The 
main one is associated with the fact that the fluid viscous dissipation can no-longer be used as a 
regularizing mechanism for the tangential component of velocity of the fluid-structure interface, as 
is the case with the no-slip condition, where the regularity of the fluid-structure interface is directly 
influenced by the fluid viscosity through the trace of the fluid velocity at the interface. As a result, 
new compactness arguments had to be used in the existence proof, which are based on Simon’s char¬ 
acterization of compactness in i 2 (0,T;i?) spaces [53], and on interpolation of the classical Sobolev 
spaces with real exponents H s (or alternatively Nikolskii spaces A^ s,p [53]). Furthermore, to deal 
with the non-zero longitudinal displacement and keep the behavior of fluid-structure interface “un¬ 
der control”, we had to consider higher-order terms in the structure model given by the bending 
rigidity of shells. The linearly elastic membrane model was not tractable. Due to the non-zero 
longitudinal displacement additional nonlinearities appear in the problem that track the geometric 
quantities such as surface measure, the interface tangent and normal, and the Jacobian of the ALE 
mapping, which are now included explicitly in the weak formulation of the problem, and cause var¬ 
ious difficulties in the analysis. This is one of the reasons why our existence result is local in time, 
i.e., it holds for the time interval (0,T) for which we can guarantee that the fluid domain will not 
degenerate in the sense that the ALE mapping A t? (f) remains injective in time as the fluid domain 
moves, and the Jacobian J v of the ALE mapping remains strictly positive, see Figure [3] 

Degeneration of the fluid domain is associated with the “contact problem” between structures, 
as shown in Figure [3 It is well known that due to the no-collision paradox associated with the 
no-slip condition [26] [27] [52], contact between two “smooth” structures is not possible in the case 
when the Jacobian J v becomes zero, corresponding to the situation in Figure [3] right. This gives 
rise to various difficulties in the numerical simulation and modeling of problems such as, e.g., heart 
valve closure, if the no-slip boundary condition is used at the fluid-structure interface. Our analysis 
presented in this paper is a first step in the direction towards studying contact between elastic 
structures in flows with slip boundary condition, which promises to shed new light on modeling of 
various biological phenomena, including the closure of human heart valves. Further research in this 
direction is necessary. 
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